Absorbing Phase Transitions with Coupling to a Static Field and a Conservation Law 
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The coupling of branching-annihilating random walks to a static field with a local conservation 
law is shown to change the scaling properties of their phase transitions to absorbing states. In 
particular, we find that directed-percolation-class transitions give rise to a new universality class 
distinct from that characterizing the depinning of the so-called linear interface model. 
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Among the many works aiming at an understanding 
of universality in out-of-equilibrium critical phenomena, 
those on transitions to absorbing states play a leading 
role because these phase transitions have no equilibrium 
counterparts and even occur in one-dimensional (ID) sys- 
tems 0,01 ■ In the simple and general case of the reaction 
and diffusion of identical particles A without site occupa- 
tion restriction, recent numerical progress (awaiting ana- 
lytical confirmation) has led to a global picture involving 
four basic universality classes Consider reactions of 
the type [mA — > (to + k)A, nA — > (n — I) A] where to, n, 
k, and I are positive integers. Outside the prominent di- 
rected percolation (DP) class, whose simplest representa- 
tives are given by single-particle reactions (to, n = 1), the 
PCPD and TCPD classes (for "pair/triplet contact pro- 
cess with diffusion" — a vocable largely used for histor- 
ical reasons) are respectively characterized by reactions 
involving two (m, n = 2) and three (to, n — 3) particles 
0. The fourth class has been considered to be defined 
by the conservation of the parity of the number of parti- 
cles, hence its usual name PC (for "parity-conserving") 
0, but it is now clear that this is not its defining fea- 
ture Nevertheless, the relevance of conservation laws 
is ascertained at equilibrium, and it remains important 
to explore it within absorbing phase transitions (APT), 
since conserved quantities abound in physical situations. 

As a matter of fact, it was argued that DP-class prob- 
lems where the order parameter is coupled to an auxiliary 
(diffusive) field with conservation generally show non-DP 
critical properties 0] . The particular case where the aux- 
iliary field is static was conjectured recently to lead to yet 
another class of APT pi |S| • (Note that the correspond- 
ing models then possess infinitely-many non-connected 
absorbing states and strong memory effects.) Numeri- 
cal simulations confirmed only partially the above ideas. 
The situation, in our view, remains unsatisfactory for 
the static case: while in ID no definite conclusions could 
be reached in 2D and 3D the numerically-estimated 
critical exponents, when self-consistent, were found to 
roughly coincide with those of the depinning transition 
of the so-called linear interface model (LIM) . A sig- 
nificant departure from DP-scaling was thus found, but 
the existence of a separate universality class remains in 



question. Meanwhile, no rigorous result is available, al- 
though a heuristic mapping on the dynamics of LIM was 
advocated in [Tl] ]. 

In this Letter, we investigate the relevance of conser- 
vation laws to APT within the more general framework 
of four "basic" universality classes briefly recalled above, 
but restricting ourselves to the coupling to a static aux- 
iliary field: now the creation of active, diffusive parti- 
cles A is conditioned to the presence of passive, static 
particles B in such a way that the total number of par- 
ticles is locally conserved. We thus study the reactions 
[mA + kB — > (m+k)A, nA — > (n-l)A+lB], and find ev- 
idence for the existence of new universality classes. We 
revisit some of the models studied in 0, H, 0, but we 
mostly deal with two-species reaction-diffusion models 
derived from those introduced in We show the ex- 
istence of a specific "C-DP" ("DP with conservation") 
class distinct from both DP and LIM in space dimen- 
sions 1 to 3, and point at a crucial difference with the 
LIM which was overlooked in We also show that a 
specific C-PCPD class exists, while no significant depar- 
ture from TCPD and PC scaling could be observed in 
the presence of the coupling to a static field with conser- 
vation. 

Extending the notation introduced in we encode 
the rules by the order of the branching and annihila- 
tion reactions using letters s, p, t (standing for singleton, 
pair, triplet) for m,n — 1,2,3, followed by integers k and 
I. For instance, reactions [2A + B -> 3A, 2A -> 2B] 
are coded as "conserved PCPD rule" C-ppl2. Our two- 
species bosonic models are updated in two parallel sub- 
steps: A particles move to one of their nearest neighbors 
(strong diffusion), then on-site reactions take place, in- 
volving the riA and ub particles present locally. 

We first present our results on the implementation of 
the C-ssll rule studied in 0. In this case, the branch- 
ing reaction A + B — > 2A is performed for each of the 
ua particles present with probability p(ng) = 1 — l/2 nB , 
while the annihilation reaction A — > B is performed every 
A particle with probability q. Note that since p(Q) = 0, 
branching is indeed conditioned to the presence of B par- 
ticles. Unlike in y|, where the total density of particles 
is used as a control parameter, here we vary q, which 
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TABLE I: Critical exponents of the C-DP class in ID, 2D, and 
3D. The corresponding values for DP and LIM are also given 
for reference. Values for the LIM class, as well as estimates 
of re for usual DP models were measured by ourselves with 
system sizes and timescales of the same order as those used 
for the C-DP class. 
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allows to start from perfectly homogeneous initial con- 
ditions (typically one A particle on each site). At large 
q, annihilation dominates and a static configuration of 
B particles is quickly reached, while at small q frequent 
branching ensures a stationary density of A particles. 

Our numerical methodology is standard: monitoring 
various order parameters (e.g. pa the density of A par- 
ticles) , we first determine the critical point q c separating 
decay to an absorbing state from sustained activity, ex- 
pecting then an algebraic law (pa ~ t~ s ) with 5 = /3/vu, 
We record also, during these runs, the mean squared local 
gradient ((V/i) 2 ) of the interface h(x,t), where h(x,t) is 
the time integral of the local activity at site x until time 
t (h — y h + 1 whenever ua > 0, with h — initially). 
First introduced in 0] for DP-class models, this ficti- 
tious interface was shown to have "anomalous scaling" 
in the form of local gradients diverging algebraically at 
the critical point: ((V/i) 2 ) ~ t 2K . (Note that it is also 
the interface conjectured to behave like in the LIM by 
Alava and Muhoz for C-DP critical points [13.) After 
this first series of runs, we estimate the decay of the 
(stationary) order parameter with the distance to the 
previously-estimated threshold, yielding an estimate of 
exponent (3. Finally, exponent z is estimated via the 
finite-size scaling of the lifetime of activity at threshold. 

Table [I] summarizes our findings for the C-ssll rule de- 
fined above in space dimensions 1 to 3. Since it is impos- 
sible, here, to present all the data which led to this table, 
we only show the most significant results, while a full 
report will appear elsewhere fl3j ]. Two conclusions can 
be drawn. Firstly, we confirm that the coupling to the 
static particles B together with the conservation law does 
break the DP class. Our estimates in 2D and 3D are con- 
sistent with those published earlier while in ID 
(Fig.^l, our data exhibit more than 3 decades of scaling, 
leading to our estimate ^d DP = 0.14(1) ^ <5f p = 0.1596 
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FIG. 1: ID data from the C-ssll rule implemented as in 
[8| and from the Leschhorn automaton for the LIM 
Near-critical decay from fully active homogeneous initial con- 
ditions (system size 2 22-26 sites), (a): density of active 
particles p vs time (top 3 curves: Leschhorn model with 
p = 0.80085,0.8008,0.80075; bottom 5 curves: C-ssll model 
with 1 - q = 0.82858, 0.8286, 0.82861, 0.82863, 0.82865). Top 
inset: p x t ' 14 vs time. Bottom inset: p x f 0125 V s time, (b): 
growth of y 2 — ((V/i) 2 ) during the same runs at criticality for 
the C-ssll model (1 — q c — 0.82861, top curve) and for the 
Leschhorn automaton (p c = 0.8008, bottom curve). 



(runs showing an effective decay exponent close to the DP 
value eventually become clearly subcritical). Secondly, 
our results also lead to ruling out LIM-class scaling: our 
estimates of exponents S, (3 and z for this C-ssll rule 
cannot fully rule out LIM values (Table HJ : no discrep- 
ancy could be registered in 2D, while in ID and 3D we 
record significant but small differences (insets of Figs.Q] 
and[3J). On the other hand, it is clear from our results on 
the "roughness exponent" k that there are fundamental 
differences between a LIM interface and the fictitious in- 
terface constructed at the C-ssll critical points: in ID, 
the LIM interface shows anomalous scaling, but k is more 
than twice smaller than for the C-ssll rule (Fig.^J)). In 
higher dimensions, LIM interfaces are not locally rough 
(k = 0), i.e. there is surface tension (Fig.[2J). Meanwhile, 
we observe, to numerical accuracy, 2k = 1 — 5 ^ for the 
C-ssll critical points, as observed for ordinary APT [lif . 

In order to assess the robustness of the above results, 
we also studied a number of other rules whose "uncou- 
pled" version have APT in the DP-class, implemented in 
the spirit of [3j. Consider for instance rule C-spl2, encod- 
ing reactions (A + B -> 2 A, 2 A — > 2B). For each of the 
[tm/2J pairs of active particles present locally, either the 
branching or the annihilating reaction takes place with 
probability p or q = 1 — p. Here p is independent of 
riB, but the branching reaction is further subjected to 
the presence of enough B particles among those present 
initially. Finally, if tla — 1, only the branching reaction 
can occur (with probability p and if ns > 0) . 

Rule C-ssll, implemented as above, was investigated in 
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FIG. 2: Same as Fig. □ but in 2D and 3D Top left: near- 
critical decay of p in 2D for C-ssll rule (top curve, 1 — 
q — 0.65298, 2 14x2 sites) and Leschhorn automaton (bottom 
curve, p = 0.64173, 2 15x2 sites); Top right: same in 3D for C- 
ssll rule (bottom curves, 1 - q = 0.57908 and 0.57909, 2 llx3 
sites) and Leschhorn automaton (top curves, p = 0.52722 and 
0.527225, 2 10x3 sites); Bottom panels: y 2 = ((Vft) 2 ) for the 
same runs (left: C-ssll rule, right: Leschhorn automaton). 




FIG. 3: Same as Fig. U,, but for C-DP rules implemented 
following 0| (see text). Only near-criticality runs are shown. 
From top to bottom: C-ssll rule (p = 0.720775, 0.720785, 2 26 
sites); C-ss21 rule (p = 0.75905, 0.75906, 2 24 sites); C-spl2 
rule (p = 0.7094, 2 22 sites); C-psll rule (p = 0.83811, 2 24 
sites). Inset: p x i° vs time for the same data. 



1, 2, and 3 dimensions, yielding exponent values equal, to 
numerical accuracy, to those reported in Table [I] Rules 
C-ss21, C-spl2, and C-psll were also studied thoroughly 
in ID with, again, exponent values fully compatible with 
those found for rules C-ssll (Fig.[3J|. All these results give 
credence to the existence of specific C-DP universality 
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FIG. 4: Same as Fig. but for C-PCPD rules, for systems 
of size 2 24 or 2 26 sites. Top 5 curves: C-ppl2 rule (from 
bottom up: p = 0.81813, 0.81814, 0.81815, 0.81816, 0.81817). 
Bottom 5 curves: C-pp22 rule (from bottom up: p = 0.82305, 
0.82308, 0.82309, 0.8231, 0.82312). Top inset 
time. Bottom inset: p x t ' 2 vs time. 
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class, distinct from both the DP and the LIM classes. 

We now report on our results on conserved PCPD, 
TCPD, and PC systems, which we only investigated in 
ID. For rules C-ppl2 and C-pp22, our estimates of critical 
exponents reveal an identical departure from the PCPD 
class which we interpret as testifying to the existence of 
a C-PCPD class. In particular, the critical decay ex- 
periments allow to clearly rule out the PCPD value of 
the 5 exponent (Fig.^J. We find, for both rules studied: 
5 = 0.17(1), 13 = 0.32(2) and z = 1.55(5), whereas our es- 
timates for the PCPD class are 5 = 0.20(1), (3 = 0.37(1) 
and z — 1.70(5). The careful study of rules C-ttl2 and C- 
tt22 did not reveal any significant departure from TCPD 
scaling for these C-TCPD rules. Similarly, we found only 
small influence of the coupling to a static field for gen- 
eralized voter (PC) rules C-mp22 and C-mp42 0|. We 
have, at present, no reason to believe that the C-PC and 
C-TCPD classes "do not exist" . 

Before our main conclusions, we come back to the sug- 
gested equivalence between the LIM and C-DP classes 
which is contradicted by our numerical results. In the 
Leschhorn automaton [l5j , used here for our estimates of 
LIM class exponents, an integer-height interface evolves 
in parallel according to: 



1 iff V 2 h L +fi>0 



(1) 



where i is a lattice index and V 2 /ii is the discrete Lapla- 
cian calculated with the nearest-neighbors of i. A key 
ingredient is the quenched bimodal random "force" fi 
(fi = ±1 with probability p and 1 — p) which is redrawn 
every time the interface advances. 

The evolution of the fictitious interfaces from which 
exponent k is estimated in C-DP models can be cast into 
an equation very similar to Q JU]- This is perhaps 
best seen when considering the Manna or "fixed energy" 
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sandpile, a much-studied model of the C-DP type which, 
unfortunately, suffers from severe corrections to scaling 
0, El 0| . There, if the number of "sand grains" 
present at lattice site i is equal to or exceeds a a threshold 
z c usually taken to be 2, then z c grains independently hop 
to a randomly-chosen nearest-neighbor. This defines ac- 
tive and passive sites, corresponding to the above A and 
B particles, and the local conservation law is obviously 
that of the sand particles. The fictitious interface thus 
simply encodes the number of times a site has discharged 
its grains since t — 0: hi — > hi + 1 when zi > z c . This 
local discharge condition can be expressed in terms of 
zf 1 and z° ut , the cumulated number of particles received 
from and distributed to the neighbors. Consider, for sim- 
plicity, the ID case with z c = 2. Clearly, z° ut = 2hi, 
whereas z™ can be written as the sum of a mean flux 
(h i+ i + hi-i) and a fluctuating part r,(/ii) [llllr^ . The 
discharge condition is then expressed 

hi+i + hi-! - 2hi-l + zi(0) + n(hi) > 

where Zi(0) is the initial number of particles on site i. In 
the general case, the fictitious interface thus obeys 

hi^hi + 1 iff V 2 /i, - 1 + 2,(0) +n(h l ) > (2) 

which is, up to the presumably irrelevant "columnar 
noise" Zi(0) — 1, formally equivalent to Eq.QJ. There ex- 
ists, however, a fundamental difference between the two 
equations: whereas fi is a delta-correlated noise, Tj is a 
sum of random variables and its amplitude diverges as 
t — > oo. As a matter of fact, since zj is bounded, both 
((V 2 /i) 2 ) and (r(/i) 2 ) must diverge in the same way. We 
have checked numerically that it is the case, and found 
moreover that ((V 2 ft-) 2 ) ~ t 2K . On the other hand, for 
the Leschhorn automaton interfaces, the amplitudes of 
both the noise and the Laplacian remain bounded, indi- 
cating that it constitutes a reasonable physical interface 
model, even though ((V/i) 2 ) diverges in ID. 

These results indicate that the connection between C- 
DP and LIM via the fictitious interfaces is only formal 
and bears no reality. Thus, the numerical values of ex- 
ponents 5, z and (3 are probably coincidentally very close 
to each other in both classes, and C-DP critical scaling 
defines a universality class distinct from that of the LIM. 

To sum up, we have found evidence that the coupling 
of usual reaction-diffusion systems to a static field with 
a local conservation law changes the scaling laws of their 
transitions to absorbing states, giving rise to bona fide 
universality classes distinct from previously known ones. 
Naturally, these results call for rigorous analytical confir- 
mation which will probably only happen when the clas- 
sification of simple APT (DP, PCPD, etc.) is itself given 
some solid theoretical basis. 

Reaching a global understanding of absorbing phase 
transitions with coupling to an auxiliary field requires, 
from our viewpoint, to study the respective importance 



of the static or diffusive character of this field and of the 
presence of a conservation law. In particular, the sole 
coupling to a static field (without conservation) is usually 
believed to be irrelevant from existing work on simple 
fermionic models such as the pair contact process pL8j . 
This should now be tested against the general framework 
put forward here. Similarly, the case of a diffusive field 
with conservation should be extended beyond the DP- 
class. Both directions are currently being investigated. 

We thank M.A. Muhoz and F. van Wijland for many 
fruitful discussions. 
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